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Abstract 

I study the leading effects of discretization errors on the low-energy part of the spectrum of 
the Hermitian Wilson-Dirac operator in infinite volume. The method generalizes that used to 
study the spectrum of the Dirac operator in the continuum, and uses partially quenched chiral 
perturbation theory for Wilson fermions. The leading-order corrections are proportional to (a 
being the lattice spacing). At this order I find that the method works only for one choice of sign 
of one of the three low-energy constants describing discretization errors. If these constants have 
the relative magnitudes expected from large Nc arguments, then the method works if the theory 
has an Aoki phase for m o? , but fails if there is a first-order transition. In the former case, the 
dependence of the gap and the spectral density on m and o? are determined. In particular, the 
gap is found to vanish more quickly as — > than in the continuum. This reduces the region 
where simulations are safe from fluctuations in the gap. 
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I. INTRODUCTION 



Simulations of lattice QCD usinsr Wilson fermions are now able to enter into the chiral 
regime Q, in large part due to recent advances in algorithms S Q. An obstacle to 
further reduction is the possible presence of arbitrarily small eigenvalues of the Wilson- 
Dirac operator, due to the breaking of chiral symmetry. Such small eigenvalues can lead to 
algorithmic instabilities and slowing down.^ This issue has been investigated numerically 
in Ref. 0], where it is found that the average spectral gap is approximately proportional 
to the quark mass, with a distribution whose width is ~ a/\/V, with a the lattice spacing 
and V the four-volume of the lattice. This allows one to estimate the parameters for which 
simulations are safe from instabilities. 

These results are encouraging, but raise further questions. In particular, what are the 
leading-order (LO) effects of discretization on the spectrum which survive in infinite volume? 
Do such effects impact the regime in which light quarks can be safely simulated? This 
paper aims to address these questions, at least in part, by determining the leading effect of 
discretization errors on the infinite volume spectral density. 

In fact, it is already known that the spectral density in infinite volume is distorted by 
discretization errors proportional to a^. Such errors lead to a non-trivial phase structure in 



the regime where m ~ d^Kq^^^^ 0, Q , in which there is competition between chiral symmetry 
breaking due to quark masses (generically denoted m) and the Wilson term. In the scenario 
with an Aoki phase |8'j, it is found that the gap vanishes throughout the Aoki phase, and thus 
for a range of quark masses, rather than for a single value (m = 0) as in the continuum 
The results obtained here quantify this further, by showing how the gap approaches zero as 
one approaches the Aoki-phase end points, and determining the distortion of the spectral 
density from its continuum form. 

The above-mentioned competition between quark mass and effects can lead to another 
scenario, in which there is no Aoki phase, but instead a first-order phase transition at a 
non-zero pion mass 0, Q]. Pions lighter than this minimal mass are not accessible. One 
does not know a priori which scenario holds for a particular choice of action. Numerical 
simulations find, however, that the first-order scenario appears to apply for simulations 



^ For a summary of possible problems, see Ref. 
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with (unimproved) Wilson fermions and the Wilson gauge action for a < 0.2 fm j9[. The 
resulting minimal pion mass is surprisingly large, roughly 500 MeV at a ~ 0.2 fm. In terms 
of the leading-order prediction, (m™™)^ = a^A^, this corresponds to a non-perturbative scale 
A ^ 700 MeV, a large but not unreasonable value. The large size of this effect provides 
further motivation for studying the discretization errors in the spectral density.^ 

I use two methods to study the spectral density and gap. Both rely on partially quenched 
(PQ) chiral perturbation theory (XPT) including the effects of discretization errors. The 
first method (which I call the "primary method" ) makes less use of partial quenching, and 
appears to be on more solid theoretical ground. It turns out, however, that it fails at leading- 
order (the order I work at here) for certain choices of the parameters ( "low-energy constants" 
or LECs) in the chiral Lagrangian. These correspond approximately to the parameters for 
which there is a first-order transition. Thus the method works predominantly if the theory 
is in the scenario with an Aoki phase. Due to this limitation, it will probably not be 
applicable to simulations with unimproved Wilson quarks and the Wilson gauge action. 
It may, however, be applicable to simulations involving improved Wilson fermions and/or 
improved gauge actions, where the nature of the phase diagram is not yet clear. 

Because of the limitations of the primary method, and as a check, I have repeated the 
calculation using a different method (which I call the "alternate method"). This method 
requires more extensive use of the properties of PQ chiral theories, properties I study using 
the replica trick. I find that the successes of the first method are reproduced, but that, in 
general, its failures are not resolved, although some avenues for further study are suggested. 

The remainder of this paper is organized as follows. In the next section I give an overview 
of the methods I use to determine the spectral density. In sec. lIIH I review PQXPT including 
discretization effects. The derivation of the primary method is given in sec. IIVI and the 
results of this method are presented in sec. Discussion and conclusions appear in sec. IVIl 
Appendix [X] gives a derivation of a key result used in sec. IIVI Appendix IHl describes the 
calculation using the alternate method. 

^ One might wonder if this large value of A is consistent with the results of Refs. P, which use Wil- 
son fermions and the Wilson gauge action, and work at small pion masses and yet see no first-order 
transition. In fact, it is consistent because the minimal pion masses (« 250, 200 and 150 MeV at 
a = 0.1, 0.08 and 0.06 fm, respectively) lie below the smallest pion masses simulated. 
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II. OVERVIEW OF METHODS 



The low-energy part of the spectrum of the Wilson-Dirac operator is related to long- 
distance physics, the classic example in the continuum being the Banks-Casher relation 
between the density of zero eigenvalues and the quark condensate 0] . The basic idea of the 
method used here is to represent long-distance physics by the effective chira 
into which the effects of discretization can be systematically introduced 



Lagrangian, 
E, U, l3- The aim 



is then to convert results from the chiral Lagrangian concerning the long-distance behavior 
of correlators of bilinear operators into a prediction for the spectral density of the lattice 
operator. 



A tool for doing so was developed for the continuum theory by Refs. 
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JJ. One 



introduces a valence quark, the condensate for which has the spectral decomposition: 

{qvqv){mv) = - [ d\-^^f^^ , (1) 
J iX + my 

where pp{X) is the spectral density of the hermitian continuum operator — zp. The key 
point is that this spectral density is independent of the valence quark mass, since it does not 
enter the quark determinant. Thus the relation can be inverted (unlike the corresponding 
expression for the sea-quark determinant): 



= -27rpp(A), (2) 

mv=— iA 



SO that one needs to calculate the valence condensate for complex valence quark mass. This 
can be done by analytically continuing the expressions obtained using PQXPT [15] for real 
masses. 

The leading-order analysis is very simple ll^. Assume that the sea-quark masses are 
small and positive, so the sea-quark condensate {qsQs) is negative. Then, for real my, the 
valence condensate has the same magnitude as that of the sea-quark, but with a sign given 
by that of my.^ Analytically continuing to complex my gives 

{qvqv) = {qsqs)\fmlr/™v ■ (3) 
The discontinuity for imaginary my gives, using eq. the Banks-Casher relation 

pMA) = -^[1 + 0(|A|/Aqcd)] . (4) 



This result is demonstrated in passing in Appendix IbI using the rephca trick. 
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The NLO correction (the term hnear in |A|) was also obtained in Ref. (reproducing and 
extending earher results of Ref. ^^^1 the calculation was extended to the e-regime in 



Ref. [ij], giving the microscopic spectral density. 

My aim in this paper is to generalize this simple continuum LO analysis to the lattice 
theory. This turns out to be more involved than one might have anticipated. Using the 
most direct generalization of eq. © leads to complications which I only partly resolve in 
Appendix El This approach I name the "alternate method" . Instead, in the main text I 
follow a different, though related, approach, which I call the "primary method". In the 
following, I first present the key equation of the alternate method since this provides useful 
background for the subsequent development of the primary method. 

To generalize eq. one cannot simply replace ]Z) with the Wilson-Dirac operator D-^, 
since the latter is not anti-hermitian, and has complex eigenvalues. Instead, the spectrum 
of interest is that of the hermitian Wilson-Dirac operator 

Qm = l^iDw + mo) = Qi, , (5) 

with rriQ the bare quark mass. To generalize eq. (Q) I note that to introduce the 75 needed 
to convert the propagator {Dw + rrio)^^ to the desired Q^, one should consider the valence 
pseudoscalar condensate. Furthermore, the extra 75 in Qm means that the analog of my in 
is a valence twisted mass. I find it convenient to introduce a degenerate pair of valence 
quarks, for which one can write a flavor non-singlet twisted mass term as ifiqvldTsQv- It is 
then straightforward to derive an analog of (Q), namely 

{qvl5r3qv)W = - J dX ^— p— , (6) 

where Pq{\) is the spectral density of Qm- As before, the use of a valence quark allows one 
to invert this relationship, yielding 



Disc[(gy75^3<?y)] 



= 227r[pQ(A) + pQ(-A)] . (7) 

fi=i\ 



Note that the untwisted part of the bare valence quark mass is to be kept at the same value 



as the mass in Qm, namely tuq. The result ((Tj) is very similar to the equation used in Ref. 17 1 
to investigate Pq(A) numerically. 

To use eq. ((Tj) one needs to determine the PQ twisted condensate for real twisted mass 
and analytically continue to complex values. For an unquenched theory, and working at 



leading-order in XPT, this would simply require minimizing the potential for real masses 
and analytically continu ing the resulting expression. This minimization has been worked 



out in Refs. 
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2l|, and analytic continuation is straightforward as the equation 
to be solved is a quartic polynomial. For a PQ theory, however, the methodology is less 
straightforward. If one uses the graded-symmetry method |l5i] , the complications are largely 
due to the presence of ghosts. The similar problem of studying the phase structure in 



quenched QCD with Wilson fermions was considered in Ref. [2^, and it was only after 
making a number of assumptions beyond those of XPT that a result could be obtained. 

Because of this complication, I decided to study eq. using the replica approach to 
PQXPT [23j. This replaces the complication of ghosts with the need to analytically continue 
in the number of valence quarks, Ny. This continuation becomes non-trivial if there is phase 
structure which depends on Ny, as might be the case here. My analysis of these issues is 
preliminary, and I have only pursued the method far enough to determine the assumptions 
necessary to reproduce the results of the primary method, and to obtain some hints for 
possible problems with that method. Because this alternate analysis is incomplete, I relegate 
it to appendix iBl 

The primary method I use makes use of the formulation given in Ref. j^. In this one 
considers the spectrum of the square of the Hermitian Wilson-Dirac operator, whose 
spectral density is denoted pi^a) (without a subscript). The eigenvalues of are both real 
and positive, with a the smallest. To get a feel for p(a), it is useful to see the form it attains 
in the continuum limit. Then = — + m^, so the spectral gap is m^, and the spectral 
density is given in terms of that of —ip by 



p^aY"""' = < ~ . (8) 

I < a < m^. 

Here I have used the results that the spectrum of —ip is symmetric under A ^ —A, and 
that exact zero modes give vanishing contributions in infinite volume. Using the Banks- 
Casher relation, eq. (j3)), one sees that p{aY"^^ has an integrable square-root singularity at 
the position of the gap 



p(«)cont^ ^^^ + 0(A^Cd) «>^- (9) 

Note the first term on the right-hand-side (which is the LO contribution in XPT) dominates 



for eigenvalues satisfying < a ^ ^qcd) since for them p(a) ^ ^qcb- is for such 
eigenvalues that the discussion in this paper applies. 

In order to discuss the renormalization of the spectral density, Ref. P introduces a 
resolvent 

R{z) = r da . (10) 

The factor of a"^ insures convergence in the ultra-violet (where p approaches its free field 
form). Thus R{z) defines an analytic function of z with a cut along the real axis beginning 
at a, across which there is a discontinuity 

—2iTip{a) 



Disc[i?] 



(11) 



Assuming a non-vanishing gap, i.e. a > 0, i? can be expanded about z = 0: 

i?(z) = f^z^Mfc, M, = - r da^. (12) 

This expansion is convergent for |2;| < a. It is useful thanks to the fact that each term in 
the sum can be written as a partially quenched expectation value jsj 

Mu= Yl «'""'(^12(^1)^23(X2) . . . P„-ln(Xn-l)P„l(0))pQ = {P"") PQ , U = 2k + Q . 

(13) 

The notation here is as follows. There are at least n flavors of valence quarks, and Pjk is 
a local bare flavor non-singlet pseudoscalar density: Pjk{x) = qjl^qk^x) . The flavors are all 
degenerate, with the same bare mass, mg, as the sea-quarks of interest. The Xi represent 
lattice sites, and are summed over, so that all insertions are at zero 4-momentum. The 
choice of flavor indices implies that the quark contractions form a single cycle, so that, upon 
insertion of the eigenvalue expansion of the quark propagator, one recovers the form of 
in eq. (fT^ . The shorthand notation {P"')pq is for later convenience. 

In order that the spectral density has a well-defined continuum limit, it must be renor- 
malized. This is accomplished straightforwardly by renormalizing the coefficients in the 
expansion of the resolvent ^j: 

oo 

Rr{z) = ^'Mk,R , Mk,R = Zf+^Mu = {P^r'+')pq . (14) 

fe=0 

Here Zp is the renormalization factor for the non-singlet pseudoscalar density (defined in a 
convenient scheme), and Pr = ZpP. It then follows that the renormalized resolvent is given 



in terms of the renormalized spectral density by an equation of the same form as eq. ()lUp. 
and has a discontinuity of the same form as eq. (fTTj) . but with p(a) — > p_r(q;), where 



Pi? (a) = Zp p{Zpa) . 



(15) 



The point of this result is that it shows how to rescale the argument and normalization of the 
lattice spectral density p(a) in order to obtain a quantity which has a good continuum limit. 
In particular, when a — 0, pnia) becomes the p(a)™'^* discussed above. This completes the 
review of the results of Ref . . 

The strategy of my primary method is straightforward. First, calculate the correlation 
functions (-Pr)pq using PQXPT for Wilson fermions (PQWXPT), and thus obtain the co- 
efficients Mk^R. Second, sum the series in eq. ()14j) within its domain of convergence. Third, 
analytically continue the result and determine the discontinuity across the real axis, thus 
obtaining the renormalized spectral density and its gap. One can relate the result back to 
the bare lattice spectral density and gap using eq. (fTKj) . 

In practice, one can work only to a given order in chiral perturbation theory. Here I work 
only at leading-order — tree-level — which is simple enough that one can sum the series. As 
I show, this order reproduces the Banks-Casher result in the continuum, i.e. it gives the 
dominant contribution for low eigenvalues. The new feature here is the extension of this 
result by the inclusion of discretization effects. 

In order to facilitate the calculation, it is useful to introduce an auxiliary function. 



which includes contributions from (-Pr)pq starting at ?7, = 2 rather than n = 6,^ and in 
which there is an additional power of i in the denominator. Given F{z), the resolvent can 
be obtained from 



^ Strictly speaking, to renormalize F{z) one needs to subtract the mixing of {Pfi)pQ and {P^)pq with the 
identity operator, mixing that is quadratically and logarithmically divergent, respectively. This mixing 
appears in XPT through contact terms between sources. In practice, one can ignore this subtlety (by 
dropping these contact terms), as the £ = 1,2 terms in the sum play no role in the development of the 
non-analyticity that is the focus here. 



oo 



.e 




(16) 




(17) 
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Thus the definition of the resolvent, eq. (fTUjl . can be viewed as a twice subtracted dispersion 
relation for F'{z). It follows from eq. ()17|) and (fTT|) that F'{z) has a cut along the real axis 
in the same position as Rr{z), but with discontinuity 



Disc[F'] 



-2t7rpR{a) . (18) 



F{z) can be obtained from F' by integrating radially outward from the origin. It follows 
that F is analytic except along the same cut on the real axis, with its discontinuity being 



Disc[F] 



pa 

= —2iTcNji{a) = —2m I da' pR{a) 

-a J a 



(19) 



as long as the integral exists, as will be the case here. The quantity N^la) is the integrated 
spectral density: the number of eigenvalues per unit volume with eigenvalue less than a. 
Thus, if one can calculate F{z) by summing the series in eq. and analytically continue 
to determine its cut, one directly obtains Nji^a). 



III. PARTIALLY QUENCHED CHIRAL PERTURBATION THEORY FOR WIL- 
SON FERMIONS 

In this section I review the ingredients needed for a calculation of the correlation functions 
{Pr)pq using PQWXPT.^ I work in the regime in which the expansion parameters p^/A^ 
(with = 47r/7r), m/AqcD (with m a renormalized valence or sea-quark mass) and a^Ag^p 
are all comparable. I assume these quantities are small compared to unity, so that it is 
reasonable to keep only the leading term in the joint chiral-continuum expansion. 

I use the graded-symmetry formulation of PQXPT, although, as will become apparent, the 
calculation is step-by-step equivalent to that using the replica trick. The partially quenched 
chiral Lagrangian including discretization effects up to O(a^) has been given in Ref. jl^ . 
The LO terms are^ 

= ^Str (9^S9,.St) - :^Str (x'+S^ + ^X-) 

- a'W^ [Str (S + T))f - a'W'^ [Str (S - St)]' - a^W'l Str {Y? + [Et]^) (20) 



^ I assume that PQXPT is a valid effective ttieory describing PQ simulations. This issue has been discussed 
in Ref. i2|. 

6 The relation to the constants of Ref. [13 is as follows: = W^-Wfi+ L^, W!j' = - W7 + Lr and 
W^' = W^-Ws+Ls. 
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Here the field S G SU{Ns + Nv\Nv) contains the pseudo-Goldstone bosons and fermions 



of the graded-symmetry group SU {Ns + Nv\Nv) l x SU {Ns + Nv\Nv) r |l5|,|2|,and "Str" 



indicates supertrace. The LECs which appear are as follows: /, the decay constant in the 
chiral limit, normalized so that = 93 MeV; Wq ~ ^qcd' which appears in the definition 
a = 2Woa, and represents the generic size of discretization effects; and the dimensionless 
constants W^'^g, which describe variations from this generic size.^ 

The quantities x± gQ- pHj) contain the renormalized quark mass matrix and the (ma- 
trix) source for the renormalized pseudoscalar density: 

x'± = 2So(M'±p), (21) 

Here Bn is the standard continuum LEG. The primes on x' and M' follow the notation of 



Ref. |25| and indicate that the 0{a) term has been absorbed by a shift of the quark mass. 
The expressions (j2H) are somewhat unconventional — usually one has x'+ = 2Bq{s + ip) 
and x'- = {x'+)\ with s and p hermitian matrix sources, the former ultimately set to 
M'. This corresponds at the quark- level in the continuum effective Lagrangian to using 
= ^(■^ + '^l^P)l) which is the natural choice of the source in an operator treatment since 
it is hermitian. What is of interest here, however, are matrix elements of Qj'^^qk without the 
factor of i. These are obtained if the source term in the continuum quark-level Lagrangian 
has the form Cq = q{M + 75p)g (with p not necessarily hermitian). A standard spurion 
analysis then leads to the result of eq. (PT|) . No problems are caused by having a complex 
action since the sources are being used as a tool to develop perturbation theory. Note that 
one obtains the desired pseudoscalar density using 

P^kA^) = = -^(S^ - S)., , (22) 

dpjk{x) 2 

where S is the action. This is the same form as in the continuum — discretization errors lead 
to corrections suppressed by powers of a. As in the continuum, one automatically obtains 
at LO the renormalized pseudoscalar density from the XPT calculation, with the scheme 
dependence implicitly contained in that of the constant Bq. 



^ For my present purposes, this notation is redundant, as Wq can clearly be absorbed into W^'_g. I use this 
notation, however, in order to be consistent with Ref. [l^ . 
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The quark mass matrix is given by 

M' = diag(m^t, m^, m^, my, my, . . ^ my, my, . . , (23) 

Nvterms Nvterms 

where I have speciahzed to A^^ = 3 hght sea-quarks. All masses are renormalized, and 
are related to the corresponding bare quark masses by an additive shift of size 1/a, an 
overall renormalization, and the further additive shift of 0(a) mentioned above j\. To 
obtain the spectral density for one of the hermitian Wilson-Dirac operators entering into 
the determinant my must equal one of the sea-quark masses. Since it is the light quark 
determinant that has the smallest gap, and since almost all simulations to date work in the 
isospin-symmetric limit ("2+1 flavors"), I choose my = m^ = m^ <^ m^. I imagine the 
strange quark mass held fixed (at or near its physical value), while the common light quark 
mass is extrapolated towards its physical value. In fact, this choice allows the strange quark 
to be integrated out of the chiral theory, avoiding concerns about the convergence of XPT 
if it is kept. In the following I assume that this has been done. Thus I will be considering 
the effective theory with 2 light degenerate sea-quarks, for which the effective Lagrangian 
has exactly the same form as eq. (j^UI) . except that (i) all matrices lose one row and column, 
(ii) ms does not appear in M', and (iii) all the LECs have an implicit, unknown dependence 
on ms- Since mg is ultimately to be fixed to its physical value, and since the constants are 
unknown, the dependence on does not add further uncertainty. Note that the subsequent 
considerations will also apply directly for a two-flavor theory (i.e. one in which the strange 
quark is quenched) — the only change being that the LECs will have different values. 

This completes the description of the PQ chiral theory. Before proceeding to the calcu- 
lation, I make some general comments on the set-up. First, the result that the LO term 
linear in a can be absorbed into the quark mass immediately implies that the discretization 
errors in the spectral density will be of O(a^). Note that this does not mean there will be 
no errors linear in a: such errors are present, but of next-to-leading-order (NLO), because 
they are multiplied by m or p^. 

Second, the LECs, and in particular IVg-S' independent of Ny. This is clear for 
those linear combinations that remain when one considers the unquenched SU{2) subgroup 

® This follows from the discussion after eq. 1)13(1 . It is of interest also to consider mv 7^ rns, with ms a 
generic sea-quark mass. This tells one about the spectral density for partially quenched quarks. I have 
not considered this possibility in this work. 
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(discussed in more detail below). The effective theory in this subgroup is independent of 
For other linear combinations, which contribute only to PQ correlation 



functions, the argument is more indirect. These combinations give the LO contribution to 
particular amputated PQ correlation functions (which can be chosen so as to cancel the 



contributions of the physical LECs, along the lines followed in Ref. j27[). At the quark- 
level, these correlation functions are independent of Ny-, by construction. Since one can, 
in principle, make LO XPT be arbitrarily accurate by sending m, — * 0, the "unphys- 
ical" combinations of LECs must also be independent of A'^v'- One can also understand 
this a posteriori: the spectrum of Qm, which knows nothing about iVy, depends on these 
combinations. 

The third general comment concerns the relative magnitude of Wq_^. It is straightforward 
to show that, for (the number of colors) large, PVg/^s ~ W^f/Ws ~ l/^c Thus one 
expects, even for A^^c = 3, that the dominant LEG describing discretization errors will be 
W^. I will make use of this expectation when describing the implications of subsequent 
results, although the results themselves will be general. 

Finally, I note that the partial quenching needed for this calculation is, in some sense, 
"minimal". For one thing, since my = ms, there are no double-pole contributions to 
propagators. While these would not appear until NLO in the present calculation, their 
absence removes the most unphysical aspect of PQXPT. In addition, in the LO calculation 
I carry out below, the ghosts (or replicas if using the replica trick) play no role. The 
information that PQXPT is "supplying" is that (a) the mass of pseudo-Goldstone bosons 
(PGBs) composed of valence quark and antiquark is the same as that of PGBs composed of 
sea-quarks [which follows from the exact SU{Ns + Ny) vector symmetry, assuming this is 
unbroken] and (b) that one can separate the contributions from individual contractions in an 
effective field theory framework. One cannot avoid partial quenching, however, because it is 
not possible with sea-quark correlators alone to pick out the quark-connected contractions 
that contribute to {Pr)pq- 

I close this section by recalling the essentials of the phase structure in the two flavor theory 
caused by the terms. ^ The phase structure is determined by minimizing the potential in 
the physical sub-sector of the theory, i.e. that involving only the two sea-quarks. Restricting 



The effects of higher order terms have been considered in Ref. |28|, |29| . 
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S to this sector, and setting the source p to zero, one finds 

VuQ = -^|^Tr(St + j:)-a'W' [Tr(St + S)]' . (24) 

Here I have used the fact that S G SU{2) to simphfy terms, defined W = + /2 
(following Ref. 2^), and dropped an irrelevant constant. 

At this stage it is useful to introduce a more compact notation. Extending the notation 



of Ref. 



I define 



m = 2Bom , wq = — — — , Ws = — — , w = — — — = wq + — . (25) 

The phase structure depends on the sign of w' [7!|. (The relation to the notation of Ref. [3| is 
w' = — /^C2.) For w' < there is an Aoki phase, running between end points at m = ±2|w'|. 
The presence of a non-zero twisted condensate inside the Aoki phase implies that the gap in 
the spectrum of Qm vanishes 0. Thus the methodology developed here apphes only away 
from the Aoki phase. In this region, the pion mass is given by 

ml = \m\ + 2w' , (26) 

which vanishes at the end points. 

For w' > there is a first-order transition at m = 0. The result ()26|1 still holds, showing 
that there is a minimum pion mass proportional to a. The twisted condensate vanishes 
along the entire Wilson axis, implying that the gap is non-zero for all m. 



IV. DERIVATION OF PRIMARY METHOD 



The aim in this section is to use PQWXPT to obtain an explicit expression for the 
auxiliary function F{z), defined in eq. (fT^ . To do so, one needs to determine {Pji)pq in a 
sufficiently explicit form that one can sum the series defining F. I will be able to do this 
at LO in the chiral expansion, with the result, after a rather extended saga, being given in 
eq. (jSOD- 

The first step is to map the lattice correlation function into one in the continuum effective 
theory. The former involves products of sums of the form These map into 

integrals J d'^x Pij^R{x), up to corrections which I will now argue are suppressed by powers 
of a, and thus are of higher order. Corrections arise from three sources: (I) in the mapping 
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of the local operators from the lattice to continuum effective theory (i.e. the usual 0(a) 
corrections to local operators in the Symanzik improvement program); (II) in the mapping 
of sums into integrals; and (III) from a mismatch between contact terms in the lattice and 
continuum theories. (I) and (II) lead to corrections to the local continuum operators of the 
form Pr{x) Pji{x)[l + 0{am)] +a^05(x) where O5 is a dimension 5 operator transforming 
as a pseudoscalar (e.g. q'-y^D'^q). Since there is no chiral suppression of the coupling of the 
pseudoscalar density to pions, and correspondingly no enhancement of the coupling of O5, 
the corrections are suppressed by the explicit factors of a. The third source is more subtle, 
but can be seen to be suppressed by 0{aml)}^ The conclusion is that, at LO, I can use the 
naive transcription 



and calculate the right-hand-side in PQWXPT. I will drop the subscript "cont", as the 
subsequent calculation is entirely in the continuum. 

From now on, without loss of generality, I take the quark mass to be positive. It must 
be large enough to avoid a possible Aoki phase (where we know that the gap closes and the 
method breaks down). The condensate in the sea sector (which I will also refer to as the 
unquenched [UQ] sector) is then aligned in the standard way, (Sc/q) = 1, and, in particular, 
parity and flavor are not spontaneously broken. Assuming the graded vector symmetries are 
unbroken this further implies that (Spg) = 1, i.e. that the partially quenched condensate 
is unity '24]. One then expands in pseudo-Goldstone bosons and fermions in the usual way, 
S = exp(2i$//). Since parity and flavor symmetries are unbroken, = for n odd. 

It is useful to look at low-order examples to understand the structure of the contributions. 
For brevity, I refer to all pseudo-Goldstone particles as pions. The diagrams contributing 
at tree-level to (-Pr)pq for n = A and 6 are shown in Figs. Hand El respectively. The cyclic 
arrangement of flavor indices in eq. ()13|) implies that there are no disconnected contractions. 

Because of the flavor indices, the operator product expansion in the continuum is scheniaticaUy 
P{x)P{0) ^ S{x)/x^. Integrating over a ceU of size gives / P{x)P{Q) ^ aS{Q). In the lattice theory 
one also has that a'^P{0)P{0) ~ aS{0), but the coefficient of a is different. Thus the difference between 
lattice and continuum contact terms is also aS{0). The matrix elements of aS(0) are suppressed relative 
to those of J P{x)P{0) by am^, with the arising from the loss of one pion propagator. Contributions 
from the regions where three or more pseudoscalars come into contact are even more suppressed. I thank 
Martin Liischer for pointing out the presence of the contributions from contact terms. 




n-l 



{PuM^i) ■ ■ ■ Pn-ln,R{Xn-l)Pnl,R{^)) PQ,cont ) 



(27) 
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(a) 



(b) 



FIG. 1: Diagrams contributing to {Pj^)pq with n = 4 at tree-level in chiral perturbation theory. Filled 
circles represent the pseudoscalar operators. 




(a) (b) (c) (d) (e) 

FIG. 2: Diagrams contributing to {Pji) pq with n = 6 at tree-level in chiral perturbation theory. 

Since all external 4-momenta vanish, the propagators are simply [with given in 

eq. (|26p]. while the vertices arise from the terms in the LO potential [all but the first term 
in eq. (plj) ]. Note that the pseudoscalar operator, given in eq. (|22I), can produce any odd 
number of particles. 

The cyclic arrangement of flavor indices impacts the calculation in two other ways. First, 
it constrains the contribution from the two-supertrace operators. In order to get a non- 
vanishing result, each vertex can only contain a single-supertrace, as illustrated by the 
examples in fig. |21 This happens automatically for the single-supertrace terms (those pro- 
portional to x± and W^), but not for those with two supertraces (those proportional to 
and W^l). For the latter, quark-level contractions such as those in fig.^J^a-b) are not allowed 
by the external flavor indices. Only one of the two supertraces in each such vertex can be 
contracted with pions; in the other supertrace, S and E"!" must be replaced by unity. This 
implies that the term does not contribute, and that the effect of can be included by 
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FIG. 3: Examples of allowed quark-flow diagrams, with filled circles representing the pseudoscalar fields, 
(a) Quark-flow diagram corresponding to flg.^. (b-d) The three types of quark-flow diagram corresponding 
to fig. 13;. Examples of flavor labels are shown; they can also be cyclically permuted (1 — > 2, 2 ^ 3, etc.). 
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(a) 



V 

(b) 



V 

(c) 



FIG. 4: Examples of quark- flow diagrams containing vertices with two-supertrace operators, in which both 
supertraces are contracted with pion flelds. Filled circles represent pseudoscalar fields and dashed boxes the 
vertices, (a) Disallowed contribution to {Pp)pQ- (b) Disallowed contribution to {Pjii pq corresponding to 
the diagram in fig. H^a). (c) Allowed one-loop contribution to {Pr)pq- 

shifting the mass term. In particular, to calculate (P^) pq one can replace with 

^aux ^ ^ Str([m-g+p]St + S[m-g-p]) ^ ^^^^ Str (S^ + [St]^) ^ ^ ^^^^ 
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where 



m 



eff 



m + 2^6 , (29) 



and p = 2Bop. The kinetic term has been excluded because all momenta vanish in the tree- 
level computation. This prescription of replacing with V^"^ works for all n, including the 
two-point pseudoscalar correlation function (i.e. n = 2). 

I note in passing that the use of the shifted mass fh'^^ does not simply generalize beyond 
tree-level. At NLO, two-supertrace operators have contributions in which both supertraces 
are connected to pions. An example for the two-point function is shown in fig. Efc). 

16 



It is useful as a check of subsequent manipulations to determine the explicit result for 
{Pr) pq LO for the simplest cases. I find 

{P'h)pq,lo = -fiBof^ (30) 

{P'k)pq.lo = -/'(-Bo)''i|^ (31) 

(P|)p«xo = -/Wi^^5^^g^. (32) 

with = nf^ + w^, consistent with eq. (j2SI)- 

In the Aoki-phase power counting, m ~ m'^^ ~ ~ rn^, and it is easy to see that 

{Pr)pq - f\B,rm'--\l + 0{m/kl)] . (33) 

The fact that the order of the infra-red divergence increases with n shows that the conver- 
gence of the sum defining F{z) [eq. (fT^ ] appears to break down when z ~ m^, corresponding 
to the appearance of non-analyticity. Note that the NLO chiral corrections from all terms in 
the sum are equally sub-dominant in this regime. For example, although the NLO correction 
to, say, (-Pr)pq is of the same absolute size as the LO contribution to (-Pr)pq, the overall 
factors of z make the LO contributions to F{z) of these two quantities comparable in the 
range of interest. The only exception to the dominance of the LO contributions is if they 
turn out not to give a non-analyticity, while the NLO terms do. This may in fact be the 
case here for certain values of the LECs. 

Further analytic progress is facilitated by observing that the PQ correlation functions 
can be expressed in terms of simpler unquenched correlators, if one works only at LO in the 
chiral expansion. Explicitly, I claim that 

2{n - l)!(P^)pQ,LO = {{Pr,i2 + Pr,2iT)conn^lo , (n > 1) , (34) 
where I use the shorthand notation 



{{PR,12 + PR,2lT)cONN,LO = ^{^j rf'a;[P^,i2 (x) + Pp,21 (x)] | ) 



CONN,LO ' 



(35) 



In words, the new correlation function involves n insertions of the same operator. This 
operator contains only two flavors, which can be chosen to be the flavors of the sea-quarks 
since the masses are the same. Thus it is an unquenched correlation function. It is to be 
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evaluated at tree- level, using the auxiliary chiral potential V^™ of eq. keeping only the 
connected contributions. In this calculation, one can replace the supertraces with traces, 
since only two flavors enter. Finally, note that I have integrated all n of the operators 
over space-time, rather than just the flrst n — 1, as in eq. (j27p . and compensated by the 
overall factor of with V the space-time volume. This makes the expression more 

symmetric, which is useful below. The use of a flnite, rather than inflnite, volume has no 
impact on the correlation functions as long as the box size is much larger than the pion 
Compton wavelength. This is because the correlation function is connected, and thus falls 
off exponentially as the separation between operators grows. Since the volume can be taken 
arbitrarily large, flnite volume effects can be made arbitrarily small. 

I give a demonstration of eq. (j34j) in appendix ^ Here I attempt to make the result 
plausible. One should keep in mind that it is a useful trick, and not a fundamental result. 
In particular, the equality does not hold either if there are terms in the Lagrangian with 
two (super)traces, or beyond tree-level. 

The result is trivially true for odd n, since both sides of flH4j) vanish. It is also trivially 
true for n = 2, as the two sides are identical, taking into account that, on the right-hand side, 
only the cross-terms (-Pr,i2-Pr,2i)lo contribute, and there are two of these. The flrst non- 
trivial case is for n = 4. At this order there are potential disconnected contributions to the 
right-hand side (schematically of the form {Pr^i2Pr,2i)\o) i which have no correspondents 
on the left-hand side, but these are to be removed by deflnition. As for the connected 
contributions, the key point is that, even though the flavor indices on the right-hand side do 
not themselves force the contractions to form a single quark-loop, the fact that the vertices 
have a single (super)trace does enforce this. For example, the quark-flow diagram for flg. ^ 
is given by that of flg. Ot, with flavor labels changed as follows: 3 — > 1 and 4^2. The lack 
of two-(super)trace vertices is crucial here: such vertices would lead to contractions such as 
that of flg. on the right-hand side of eq. (j34|l . but not on the left-hand side. The counting 
of contractions is, of course, different in the two expectation values. On the right-hand 
side of fl34j) . one has an extra choice not available on the left-hand side (whether to start 
with Pr,i2 or Pr_2i in a particular term), leading to the overall factor of 2. In addition, the 
subsequent choices of which P to contract with are flxed by the flavor indices on the left- 
hand side, there is no such constraint on the right-hand side, so that there are (for = 4) 
an additional 3! contractions. This is most clear for flg. Q^, but also holds for flg. ^p. 
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This argument generalizes straightforwardly for the simplest diagram at each order, that 
with a single vertex producing n pions (e.g. figs. ^ and|2K): the relative contraction factor 
is 2 (n — 1)!. I show in appendix 1X1 that this result holds for all diagrams. 

Another check on the relative factor can be obtained as follows. Assuming that the two 
sides of eq. are proportional (which is established in the appendix), and that one is 
working at small enough so that the LO result is accurate, one can calculate the relative 
size of the two sides at the quark-level. Since the PQ correlator involves only a single, cyclic 
contraction, it must be that taking the connected part of the unquenched correlator (which, 
recall, is to be done at the pion-level) is equivalent to taking the connected part at the quark- 
level (again, when working at LO in XPT). This is because the quark-connected contractions 
are identical on the two sides, while quark-disconnected contractions contributing to the 
unquenched correlator have different dependences on the volume. However, the number of 
quark-connected contractions for the unquenched correlator is 2{n — 1)!, compared to only 
one for the PQ correlator, and so the relative factor is as claimed. 

To proceed it is useful to introduce a further auxiliary function: 

= E —^(iPR,i2 + PR,2ir)coNN . (36) 

n=l 

Both even and odd powers of n are included in the sum, even though the latter vanish. This 
implies that G is really a function of rather than fi. The connected correlation functions 
in the sum are to be evaluated with the auxiliary Lagrangian rather than C^. Note 
that they are not truncated at LO in the definition of G, although in practice I will work 
only at LO. For this reason I do not need to specify counterterms. 

Given the result (jH^ the functions G and F are related at LO as follows: 

^('")^« = E ^^2(n - 1)! {PJ^)PQ,L0 (37) 



n=l 



oo 



= E ^ (Pn)pQ,LO (38) 

n=l ' 

oo / 2\l 

= E {Pr)pq,lo (39) 

1=1 ^ 

= F(V)lo- (40) 

On the penultimate line I have used the fact that terms with odd powers of /i vanish. Thus if 
one can sum the series defining G{iJ,) at LO in the chiral expansion, then, by continuing from 
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positive to negative /i^ (i.e. from real to imaginary /i) one can obtain the desired function 
F in the regime where its discontinuity gives the integrated spectral density. 

Summation of the series defining G{^) is straightforward, using the fact that the log of 
the partition function generates connected correlation functions: 

G{^^) = V ^-^^{{Pr,i2 + Pr,2iT)conn (41) 



n 

n=l 



oo 



^ n\ V 

n=l 



6 6 



hl2 Sp2l' 



n 



lnZ^"^(p) 



(42) 



p=0 



= ^in[^r(/^)/^r(o)] ' (43) 

where Z'^^{^) is the partition function evaluated using the auxiliary chiral Lagrangian but 
now with the addition of a twisted mass 

^r(/^) = ^^^(a^Ea^St) + V^'^^(/i) (44) 
V-^(/i) = -f I ^i^^ ^ ^ + \ ■ (45) 

Here /I = IB^^i. Again, counterterms are not shown since I will evaluate this partition 
function only at LO. Note that the twist is in the T\ direction, because of the choice of flavor 
indices in eq. (jHlI). 

The final step uses the result that the LO partition function is given by the classical 
result: 

^r(/^)^o = exp \-S^M\ , (46) 
assuming that is real so that the action, 5*, is real. Thus one obtains 

F{-^^)lo = G{fi)Lo (47) 
= :^ln[Z-^(/x)/Z-^(0)] (48) 
= -V;^in(/^) + VZ(0). (49) 

The second term is /i-independent constant which does not contribute to the discontinuity 
and will henceforth be dropped. 

In summary, the method I arrive at is as follows. First, minimize the auxiliary un- 
quenched potential V^^^{fi) of eq. for real /i. Due to parity, the result in fact depends 



20 



on /i^. Second, analytically continue the result to complex z = —fi'^. Finally, determine the 
discontinuity for positive real z (negative /i^), which, using eq. (fTI?j). gives the integrated 
spectral density and along with this the gap: 



Disc [V^^^ifi' = -z)] 



-2iTiNR{a) . (50) 



I show in appendix ini how, under certain assumptions, one obtains an equivalent result using 
the alternate method based on eq. 

V. RESULTS 

In order to test the method, I first check that it reproduces the LO continuum result. 
Setting = = 0, the auxiliary potential is identical to that of the unquenched continuum 
theory with a twisted mass, i.e. 

This is minimized by 

E„ = 4±*t, (52) 
so that the minimum of the potential is 



V'^'^^li)^^, = -fWfn'+l^' = -2fBoVm^ + i^^. (53) 

As claimed above, this is a function of /i^ which can be analytically continued to give, up to 
an irrelevant constant, 



F{z) = 2fBoVm^ - z . (54) 
This function has the expected analytic form, i.e. it is analytic except for a cut along the 



real axis, starting at the gap a = rn?, with a discontinuity —AipB^^/ a — rn?. Using eq. (fTTHl 
one thus finds that 

Nr^o) = '^^^J y/a-n? =^ pR^a) = — ^^"^ ^ , (55) 



7rv« — m 



which agrees with the Banks-Casher result, eq. 0, since (qq) = —BqP at LO in XPT. One 
can also expand F in powers of z and, using eq. (fT^. reproduce the results for {Pr)pq in 
eqs. dnniEl for Wq = 0. 
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Now I apply the method for a 7^ 0. Mathematically, the auxiliary potential is identical 



to that which arises when studying the phase structure of twisted mass QCD, a. 



2G 



though 



2l|. In 



- - 1) ■ m 



the coefficients are different. This problem has been studied in Refs. jisl . \li 
particular, it is shown in Ref. 2(| that the minimum of the potential (for real fi) occurs 
when the condensate is aligned in the same direction as the twist in the mass term. In the 
present case this means that S = exp{iujTi) and the auxiliary potential becomes 

— — — = —fh'^^ COS u — 'jl sin Lu — — cos{2lu) (56) 

Extremizing leads to the quartic 

wy + 2wsm'^c^ + [(m'^'^)2 + Jl^ - wl]c^ - 2wsm'^c = (m'^^)^ , (57) 

where c = coscj, and I recall that a "hat" implies multiplication by 2Bq. The extrema of 
the potential are given in terms of the roots, which I label c, by 

P rrf^ + w^c 4 

Note that, as claimed above, both these extrema and the roots themselves depend on /x^, 
despite the fact that the potential itself is linear in /i. 

In order to interpret the subsequent results, it is useful to study the physical meaning of 
the parameters than can be varied. Recall that the original PQ chiral Lagrangian depended 
on two parameters describing discretization errors {w^ and ws), the quark mass (conveniently 
packaged into m = 2BQm), and the decay constant /. The latter enters only as an overall 
factor and does not play a role in determining the size of the gap or the shape of the spectral 
density. Different linear combinations of the three other parameters enter into the physics of 
the unquenched (sea-quark) theory and the spectral density. The phase structure and pion 
properties in the unquenched theory are determined hj w' = We + Ws/2 and m, as explained 
at the end of sec. IIIH while the spectral density depends on ws and in^^ = m + 2wq. The 
situation is potentially confusing because to study the spectral density one uses the auxiliary 
Lagrangian, and this is identical in structure to that describing the unquenched theory, but 
with different LECs. 

I try here to provide a guide to the parameter space. The dependence of properties of 
the unquenched theory on the parameters is shown in fig. The method I use breaks down 
in the Aoki phase {2w' < —m). If — m < 2w' < one is in the Aoki-phase scenario (though 
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FIG. 5: Dependence of phase structure and spectral density on LECs describing discretization errors, for 
fixed m, and at LO in XPT. The Aoki phase in the unquenched (sea-quark) theory occurs in the hatched 
[red] region. Lines of constant are paraUel to its boundary, with the dot-dashed hne corresponding 
to the continuum resuh — fa. This Hne also marks the boundary between Aoki-phase and first-order 
scenarios for the unquenched theory. For a given gauge and fermion action, the continuum limit is taken 
by approaching the origin at fixed 2we/'Ws, as in the example shown by the thick arrow. The method used 
here is inapplicable in the Aoki phase, and fails for ws > (light [yellow] shading). Thus it works only in 
the unshaded triangular region. 

away from the Aoki phase itself), while ii < w' one is in the first-order scenario. If one 
takes the continuum limit using a particular choice of gauge and fermion action, and holding 
m fixed, then 2wq and Ws will both decrease as a^, but their ratio res = '^wq/ws will be fixed 
up to logarithmic corrections. Thus one approaches the origin in the figure along a line of 
slope res (so that one remains in one or other of the scenarios). The issue is what happens 
to the gap and spectral density as one does so. The actual value of the slope depends on 
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the actions, but one expects from large Nc arguments that \rQs\ < 1, as in the example in 
the figure. 

It turns out that the primary method developed in this paper only yields a prediction 
if Ws < 0, as will be explained below. One way to phrase this result is that the method 
works only if theory described by the auxiliary Lagrangian is in its Aoki-phase scenario. 
Conversely, the method fails if the auxiliary theory is in its first-order scenario. Combining 
the constraints, one finds, as shown in fig. |31 that the method only works in a triangular 
wedge of parameter space. This allowed wedge contains regions of both the Aoki-phase and 
first-order scenarios of the original, unquenched theory. If one assumes that the LECs take 
their most likely values such that \rQs\ < 1, then the allowed wedge contains only the Aoki- 
phase scenario of the unquenched theory. In other words, if \rQs\ < 1, the scenarios in the 
auxiliary theory and the original theory match, and if the former must be in the Aoki-phase 
scenario for the primary method to work, so must the latter. 

The roots of eq. (j57jl depend on two dimensionless ratios, which I choose to be 

r^^= ^ and ^ 1^ = , (59) 

ml m + 2wq + Ws m^ m% 

where ml is the value for the unquenched sea-quark pion (which is unaffected by the twisted 
valence mass). Roughly speaking, r characterizes the relative size of the contributions to 
ml due to discretization errors and the quark mass. Its varies in magnitude from in the 
continuum limit to oo at the Aoki phase end points. Its sign matches that of Ws- Thus 
for Ws < 0, its range is — oo < r < 0, while for ws > 0, the range is < r < +oo. Some 
lines of constant r are shown in fig. El The importance of these lines is that the distortion 
of the spectral density and the gap by discretization effects is constant along them. When 
taking the continuum limit, one crosses these lines, and the distortion changes. The rate of 
this change depends on the slope of the approach to the continuum (i.e. on res), but the 
quahtative nature of the approach to the continuum is independent of the slope. What does 
matter, however, is whether one approaches from positive or negative Ws- 

I first describe the results for < 0, and then explain the failure of the method for 
Ws > 0- I will show results for r = —10, —1, —1/9 and —1/99, which span the "allowed" 
region of parameter space, as shown in fig. IHl (r = —1/99 lies slightly to the left of the 2we 
axis). This sequence of r values corresponds to varying the discretization errors from large 
to very small. 
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FIG. 6: Parameter space (as in fig. El showing lines of constant r, along which the form of discretization 
errors in the spectral density does not change. Note that the line with r = 1 also corresponds to rh"^ = 0. 

It may be useful to give another characterization of the size of discretization errors. This 
is the ratio, r', of the quark mass to the half-length of the phase boundary in the unquenched 
theory, L = \w'\/Bo. Note that the boundary has the same length in both scenarios, and is 
a concrete measure of the size of discretization errors. This new measure is related to r by 

L \2w'\ |r(l + r68)| ' ^ ^ 

If vqs = 0, then r' is determined by r, and the choices of r listed above correspond to 
r' = 1.1, 2, 10 and 100, respectively. In words, these values correspond to a quark mass 
very close to the Aoki-phase end point, one half-length away, ten half-lengths away and one 
hundred half-lengths away. If rgs is non-zero, the values of r' corresponding to the standard 
set of r change, but the qualitative feature of a transition from very large to very small 
discretization errors does not. 
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FIG. 7: Real and imaginary parts of roots c for r = —1, plotted against z'. The variable z' is chosen so 
that the continuum gap occurs at z' = 1. The thickest solid curve (blue) shows the root of interest, which 
has c = 1 at z' = 0. It joins with the thinner solid curve (red) at the start of the gap, and for large z' these 
form a complex conjugate pair. The other roots are shown as dashed and dot-dashed. 

I begin by showing results for r = — 1, the value for which discretization and quark mass 
contributions are comparable. The real and imaginary parts of the roots c are shown in 
fig- 13 Since z' oc — yU^, negative z' corresponds to real twisted mass where the potential is to 
be minimized. This picks the root with c = 1 aX z' = 0, i.e. that for which the condensate is 
S = 1 when the twisted mass vanishes. For 2;' < 0, c < 1, corresponding to the condensate 
being "twisted" by the twisted mass. Analytically continuing to imaginary /i corresponds 
to staying on the same root for positive z' . One finds that this root joins with another 
and becomes two complex roots at z' ~ 0.2. This leads to exactly the expected analytic 
structure, namely a branch cut on the positive real z' axis, as one can verify by expanding 
the function in the vicinity of the position where the roots join. One can thus read off 
the gap: a ~ 0.2m^/(2Bo)^. This is significantly smaller than the result one would obtain 
ignoring discretization errors, cicont = m'^ = m^/(2_Bo)^ (corresponding to z' = 1). 

As one moves away from the continuum limit, the discretization errors in the gap become 
more pronounced. Figure p(a)| shows the roots for r = —10, for which the gap is ~ 2% 
of the continuum result. By contrast, fig. |8(b)| (r = —1/9) shows the gap moving towards 
its continuum value as one reduces the discretization errors. Very close to the continuum 
(r = —1/99), the gap is almost unity, and the roots track those for a = except close to 
the gap, as illustrated in fig. IHl 
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FIG. 8: Real parts of roots for (a) r = —10 and (b) r = —1/9. Notatfon is as in fig.EI The vertical scale 
is chosen to magnify the region of interest, so that not all roots appear in (a). 
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FIG. 9: Real and imaginary parts of roots for r = —1/99 (with notation as in fig.|3 except not all roots 
appear given the vertical scales) compared to those in the continuum limit, r — (shown as a thin solid 
[black] line). Note that the horizontal scales differ in the two plots. 

I next show how the gap depends on the pion mass-squared at fixed a. In order to 
facilitate comparison with the results of Ref. ^j, I plot the gap for Qm (which is called 
fi in Ref. j^, but which I call /iq here to avoid confusion with the twisted mass). More 
precisely, I plot, in fig. |10(a)1 2BQjjQ/\ws\ = 2BQ\Ja/w1 versus m^/\w^\. Dividing by \w^\ 
gives dimensionless quantities. If Irgsl <^ 1, the a;-axis gives the pion mass-squared in units 
of the Aoki-phase half-length. The quantities are such that, in the continuum limit, the 
result is a straight line with unit slope. 
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FIG. 10: Results for the spectral gap of Q™, /ig, plotted against rn^/\wg,\. (a) Gap in dimensionless 
units, 2_BoMq/I''^8| (thick [red] curve) and result in continuum limit (thin [blue] curve); (b) Ratio of gap to 
mJ/(2i3o), which is unity in continuum limit. Results are for < 0. 



The predicted gap is smaller than the continuum expectation. To bring this out, I show 
in fig. |10(b)| the ratio of the gap to the continuum result, / {2Bq) . This approaches unity 
for m^/jwgl ^ 1, but the discretization errors reduce it significantly as one approaches the 
Aoki-phase end point. The analytic form of the result as — is 



2Bo^iQ _ 2BoV^ 



mt 



rat 



3/2 



mt 



(61) 



Note that the vanishing of the gap at the end point does not imply that the ratio of fig. 10(b) 
must vanish too. It would be sufficient if the ratio tended to a finite constant as 0. 
Thus discretization errors lead to an additional reduction in the gap, which makes simula- 
tions less "safe" in the sense of Ref. P]. The square-root behavior is related to the mean- field 
exponents that apply at LO. In fact, if one approaches too close to the end point, higher 
order terms in XPT become important, and the exponents are corrected 0, 3fll |. 

The method also yields the discretization errors in the integrated spectral density. I show 
the impact of these errors for the standard values of r in fig. ^2 The units of Nii{a) are 
such that the continuum result (|^ is^^ 



TT 



X Nnia) = 2A/a/mcfr - 1. 



(62) 
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The analytic form as r — > — oo asymptotes, for small a/m^fj, to {'nNR)/{BQpmcfi) = (3V3/2)(v/2)^/^. 
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1 2 3 4 5 miff 

FIG. 11: Results for integrated spectral density NR{a). I plot [tt/ {Bof'^nicff)] x Nji{a) versus a/m^g, 
where rrics = m + wq/Bq = fff^ / {2Bq) . Moving from top to bottom, results are for r = — 10 (magenta), 
— 1 (blue), —1/9 (red), —1/99 (green) and r — Q (black). The latter is the continuum result if = 0. Note 
that the upper two curves do have a gap (i.e. Nji{a) only becomes non-zero starting at a non-zero value of 
a), but this is not visible in the figure. 

The figure shows what happens if m is held fixed and the lattice spacing reduced. If It^el <^ 
1^8 1; then far above the gap the integrated density is unaffected by discretization errors, 
while there are significant downward shifts at the low end of the spectrum. In other words, 
aside from the smallest eigenvalues, say those with a < 3m^, the spectrum is that of the 
continuum using the mass m. Note that this is not the quark mass which vanishes at the 
Aoki-phase end point (where m^r ^0). It is the so-called PCAC mass, which vanishes 
instead in the center of the Aoki phase. The effect of a non- vanishing wq is essentially to 
shift the spectrum clS 8b whole by an amount of O(a^). 

I now explain the failure of the method for > 0. It turns out that the nature of the 
failure depends on the sign of m'^^ = m + 2wq. As can be seen from fig. IHl the boundary 
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FIG. 12: Real part of roots for r — 1/3, for two choices of scale of z'. The root of interest is shown by the 
thick (blue) solid line. 

between positive and negative m^^ is also the line along which r = 1. 

I first consider •m'^^ > which, given that ws > 0, implies that < r < 1. For Irgsl < 1 
this corresponds to the first-order scenario. I show the real part of the roots for a represen- 
tative choice, r = 1/3, in fig. ^| This corresponds to r' = 2 if rgs = 0, i.e. discretization 
errors and the quark mass contribute about equally to the pion mass. The key point to note 
is that the root of interest, which passes through c = 1 at 2;' = 0, is real for all real choices 
of z'. This remains true throughout the range < r < 1. 

This result implies that there is no discontinuity in the potential, and, using eq. ()5()j] . 
that the integrated spectral density vanishes at LO. This is a peculiar result, particularly 
because, by sending wq, Ws (implying r —>■ 0+) one can approach arbitrarily close to 
the continuum, where the method works and one obtains a non-vanishing spectral density. 
Mathematically, what happens is as follows. Although the root of interest does not join with 
others for real z', it does so at a pair of complex conjugate values. Thus there are two cuts 
in the complex plane, as illustrated in fig. El As r — 0+, the two cuts come together at 
z' = 1, leading to a single cut on the real axis, given by the continuum result eq. As r 
varies between and 1, the ends of the cuts follow the paths in the complex plane sketched 
in fig. El They come together again when r = 1, leading to a cut along the negative z' real 
axis starting at z' = —l}"^ 

The appearance of a cut on the negative real axis for r = 1 can be understood if wq = 0. Then ?' — ^ 1~ 
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FIG. 13: Analytic structure of V^^^{z') for ws > and fh'^^ > 0, corresponding to < r < 1. The cuts 
(thick [red] Unes) are shown for the example of r = 0.1. As usual, the position of the cuts is a convention, 
while their end points are fixed. The trajectory of the end points as r changes is shown by the (blue) dots. 
The values of r chosen are 0, 0.01, 0.1, 1/3, 1/2, 2/3, 0.9, 0.99 and 1.0. Note the symmetry under r ^ ^ — r, 
Re[z'] ^ -Re[z']. 

The implication of these results is that for ws > and fh'^^ > the analytic structure 
of the root of interest, which is inherited by the potential V^^^, differs from that of F{z). 
Recall that the latter has, by construction, a cut only along the positive real axis, starting 
at a > 0. Thus the result derived above relating F to V^^^, eq. (jl^ . fails. The conclusion 
I draw is that the LO approximation is inconsistent for F{z) if Ws > 0, and higher order 
terms are needed, a possibility already discussed following eq.(jnSI)- An alternative is that 
the gap vanishes, so that the primary method is inapplicable, but this is implausible as it 

corresponds to — > +oo, which means, since m is held fixed, that one is approaching the first-order 
transition boundary. But in this case there is a non-analyticity, starting at p, = ws — Tn^, namely the end 
point of the first-order transition line. I do not know of a physical explanation of the presence of this cut 
if We ^ 0. 
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FIG. 14: Real and imaginary parts of roots for r — 2. The root of interest is that shown by the thick 
(blue) line. 

does not match onto the continuum result as a — > 0. 

The method fails for essentially the same reason when ws > and rh^^ < 0. This 
corresponds to 1 < r < oo, and, in fig. El is the region below the horizontal r = 1 line. The 
generic behavior of the roots is illustrated by the example of r = 2, shown in fig. EI Because 
in'^^ < the root which minimizes the potential for real /i is that for which c = —1 a.t z' = 0. 
For 1 < r < oo, this root does not connect with any others for real z'. The situation is very 
similar to that for •m'^^ > 0, namely there are two cuts in the complex plane. As before, 
since this analytic structure differs from that of F{z), it must be that the approximations 
leading to the relation break down. 

Despite the failure of the method for ws > 0, there is no reason to doubt the validity of 
results obtained for wg < 0. These results provide a self-consistent solution for F{z) and 
thus for the resolvent Rr{z). What I mean by this is that the series for Rfi{z), evaluated 
term by term at LO, and summed, gives an analytic function with a cut only on the real z 
axis, and no other singularities. This follows from the properties of the roots of the quartic, 
but, to be sure, I have checked this result numerically by verifying the dispersion relation 
eq. (fTUjl for real z < d. I stress that the self-consistency of the LO result does not imply 
that there are no higher order corrections to pR^a) and Rr{z). 
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VI. DISCUSSION AND OUTLOOK 



In this paper I have attempted to generahze the methods developed in the continuum 
for calculating the spectral density of —ip to the operator of interest for the lattice theory, 
Qm- The generalization necessary for the inclusion of discretization errors turned out, at 
least in the approaches I have followed, to be rather involved and only partly successful at 
leading-order. 

Nevertheless, the success of the method for < leads to a result which is relevant 
for numerical simulations, namely the reduction in the spectral gap of compared to the 
continuum expectation. Since this likely applies only for the scenario in which discretization 
errors lead to an Aoki phase, it is probably not directly applicable to the numerical results 
for the gap obtained in Ref. P]. These are obtained with the Wilson gauge and fermion 
actions, for which one is apparently in the first-order scenario. 

The result may be relevant, however, for simulations that are presently underway using 
improved Wilson fermions and/or improved gauge actions, for which the phase structure has 
not been established. With this in mind, I briefly discuss the way in which the discretization 
errors in the gap might best be seen. Looking at fig. HUl one sees that, for m^/lwgl ^ 2 
a linear fit of /ig to (or, at this order, equivalently to m) would work reasonably well, 
but there would be a non-zero, negative intercept, and the slope would be smaller than the 
continuum expectation. Significant curvature should appear below this point. If one were in 
the Aoki-phase scenario with ~ and [wgl having a numerical value close to that which 
explains the first-order transition observed with Wilson gauge and fermion actions, i.e. if 

= -a^A.^ with A ^ TOOMeV, then this "critical mass" is = 250 and 150 MeV at 
a = 0.1 and 0.06 fm, respectively. Of course, if A were reduced by a factor of 2, as might 
be achieved by using alternate actions, these values for would be reduced by a factor of 
four. 

It is important to note that, since the discretization effects discussed here are of O(a^), 
there is no a priori reason to expect them to be smaller for improved quark actions. One 
must simply investigate the issue numerically. 

On the theoretical side, the failure of the primary method at LO for > is puzzling. 
One can imagine working at arbitrarily small a and m where one would expect the corrections 
to the LO form to become arbitrarily small, and yet, apparently, they must be able to lead 
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to a significant change in the analytic structure of the resolvent Rii{a). In an attempt to 
understand this puzzle I have repeated the calculation using a second method, and, although 
some new features appear, they do not appear to provide a solution for general values of the 
low-energy constants (and in particular for Wq ^ Wi ^ 0) . 

Despite this puzzle, I stress again that the result for wg < provides a self-consistent 
solution for R^i^a), and I see no reason to doubt this result. 

Beyond the task of resolving the above-mentioned puzzle, there are a number of directions 
in which it would be interesting to pursue this work. It could be applied for non-zero twisted 
sea-quark mass, although this is perhaps more of theoretical than practical interest as the 
gap is known to be bounded by the twisted mass. Of more practical interest would be to 
determine the impact of working at finite volume. I expect the corrections to the results for 
the spectral density to be exponentially suppressed until one enters the e-regime. As noted 
in the introduction, however, fluctuations in the gap have been observed to fall as inverse 
powers of V, and it would be interesting to see whether the effective theory can reproduce 
this behavior. 
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APPENDIX A: DERIVATION OF KEY RESULT 

In this appendix I demonstrate the result in eq. (jSU, which I repeat here for clarity: 

2{n - l)\{P^)pQ,LO = {{Pr,12 + PR,2ir)cONN,LO , {u > 1) . (Al) 

As noted in the text, this result is trivial if n is odd, so I assume that n is even in the following. 
Recall that both correlators are to be calculated using the auxiliary chiral potential V^"^ of 
eq. (j^Hj) . in which there are only single-supertrace operators. 
The demonstration is based on the following results. 
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(I) Only charged pion propagators appear in {{Pr,i2 + Pr,2iT)conn,lo and {P^)pq,lo- 

Charged here means being flavor off-diagonal, e.g. "12" or "35" , while neutral means 
flavor diagonal: "11", "22" .... The argument goes as follows. Tree-level diagrams 
can always be divided into two pieces by cutting any propagator. Assume that there 
is a neutral pion propagator and cut it. Since the operators Pr are all charged, there 
must be an even number in each of the two resulting pieces to add up to an overall 
neutral combination. But an even number of Pr's can, together, produce only an even 
number of pions, which, through interactions, can only change to an even number of 
pions. Thus one reaches a contradiction. 

In fact, for (Pr)pq,lo the argument is more simple: the only set of pseudoscalar 
operators which are neutral is the whole set, but that leaves no operators for the other 
piece of the diagram. 

(II) For each pionic diagram contributing to the correlators of interest, one can unambigu- 
ously associate a finite number of planar quark-flow diagrams. 

This holds true as long as one does not have neutral pion propagators, as is the case 
here given observation (I).^'^ To obtain the quark- flow diagram one simply follows 
the flavor indices, connecting them into a single "trace" at each vertex and at each 
pseudoscalar operator. The flows are directed, with arrows indicating the flow of 
fermion number. Note that I do not distinguish between diagrams in which the flavor 
labels are permuted: these are all lumped into a single quark-flow "skeleton" , with the 
counting of possible attachments of flavor labels to be done below. 

For each pion-level diagram there are, in general, multiple quark-flow diagrams, corre- 
sponding to different orderings of the contractions between the Pr. For example, the 
pion-level diagram of flg. |2t leads to the three quark-flow diagrams of flg. [Hb-d. I have 
drawn the diagrams in flg. |21to be planar, in the sense that no quark-lines cross. This 

The potential problem with neutral propagators is that they can have "hidden" quark-disconnected con- 
tributions due to the constraint Str($) = 0. Thus one does not know how to assign the quark-flow without 
further considerations. In fact, in the present case, with all valence quarks degenerate, this restriction can 
be shown to be unnecessary, but it simplifies the argument. 

Although it is not needed for the argument, I note that quark-flow can also be traced if there are vertices 
with multiple supertraces. 
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is their canonical form, for then the quark-hne runs around the outside of the diagram, 
and the connection of flavor indices is transparent. That it is always possible to draw 
the quark-flow diagrams in this way is property of tree-level diagrams, following from 
the fact, noted above, that tree diagrams fall into two pieces if one cuts any propaga- 
tor. Thus one can cut all propagators, make each vertex and insertion of Pr planar, 
and then glue the propagators back together starting from an arbitrary point in the 
diagram, resulting in an overall planar diagram. The key point is that the vertices 
and insertions can be treated independently. 

(III) All contributions to (P^) pq.lo contain a single quark-line in the corresponding quark- 
flow diagrams. 

For each quark-line in the diagram, all of which are necessarily closed, the total 
"charge" of the Pr attached to it must vanish. Given the charges in (-P^)pq (see 
eq. W^ . this is only possible if all the Pr are attached to a quark-line, so there can 
only be one such line. Examples are shown in fig. 01 This result implies the (rather 
obvious) fact that there are no disconnected contributions to {Pr)pq,lo- 

(IV) All contributions to ((Pr,i2 + Pr,2i)"')conn,lo contain a single quark-line in the cor- 
responding quark-flow diagrams. 

This result relies on the facts that V^"^ contains only single-supertrace operators and 
that there are no neutral propagators [observation (I)]. These facts imply that there 
is no way for a quark-line to "close off" within a vertex or a propagator. (Examples 
of how this is possible with two-supertrace operators are given in fig. |3^-b.) Thus, 
if there is more than one quark-line, the pion diagram must fall into as many pion- 
disconnected pieces as there are quark-lines. This cannot happen, by definition, if the 
diagram is connected at the pion-level. 

Note that this result fails beyond LO. A one-loop example of a contribution to ((-Pr,i2 + 
Pr,2iY)conn which involves two quark-lines, and is thus absent in {Pr)pq, is given 
in fig. HSI 

(V) The same planar quark-line diagrams contribute to ((Pr,i2 + Pr,2i)^)conn,lo md 
{Pr)pq,lo- 
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FIG. 15: Example of quark-flow for a one-loop diagram involving two single-supertrace operators which 
cannot contribute to {P^i) pq but which does contribute to {{Pr.u + Pr,2iY)conn- 

Any pion-level diagram contributing to {(Pr^u + Pr,2iY)conn,lo can also contribute 
to {Pr)pq,lo, and vice- versa, since the vertices and insertions of Pr produce the same 
number of pions in both cases. Thus the pion-level diagrams are the same for both 
correlators. 

The remaining question is whether the quark-flow diagrams that contribute to a given 
pion-level diagram are the same for both correlators. That this is the case can be 
seen as follows. From (III) and (IV), all quark-line diagrams contributing to both 
correlators involve a single quark-line. Adding back in the flavor indices, these will 
alternate as 1212 ... 2 in the unquenched correlator, and as 1234 . . . (2n) in the PQ 
correlator, as illustrated in fig. ^1 Simply by interchanging these flavor labels one 
interchanges legitimate quark-line diagrams in the two correlators. Thus if a quark- 
flow diagram is present in one correlator, it will be present in the other. 

(VI) The relative contribution from each planar quark- flow diagram to ((Pr,i2 + 
Pr,2iY)conn,lo {Pr)pq,lo is the relative contraction factor 2{n — 1)!. 

The contributions to both correlators have now been broken down to those from a sum 
over (the same) planar quark-line diagrams. These diagrams are nothing more than an 
explicit representation of all the ways the flavor indices are connected when the pion 
fields are contracted. Thus all flavor factors are accounted for explicitly. The other 
factors that contribute to the diagram are the propagators, the factors which result 
from expanding out S and S"!" to obtain the vertices and insertions, and the counting 
of contractions. Only the latter differ between PQ and unquenched correlators. The 
propagators are the same because all quark masses are the same, and the expansion of 
S*^^") leads to the same factors in both theories. Furthermore, the "internal" counting 
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FIG. 16: Examples of flavor labels for a planar quark-flow diagram contributing to flg. Et for (a) {{Pr. 12 + 

PR,2lf)cONN. and (b) (P|)PQ. 



of contractions — how many ways the pion fields in the vertices can be contracted with 
each other to yield the given planar quark-fiow diagram — are also common. The only 
part of the calculation which differs between the theories is the choice of contractions 
with the external Pr. In the unquenched correlator, there are 2 x n! ways of associating 
the n operators (-Pr,i2 + -Pr,2i) with the external insertions: n\ since each operator can 
be attached to each insertion point, and 2 because, for each such attachment, one can 
globally interchange Pr,i2 and Pr,21- By contrast there are only n contractions in the 
PQ case, because cyclicity must be preserved. Thus the relative factor is 2{n — 1)!, as 
claimed. 

In summary, by breaking the diagrams up into planar quark-fiow parts, which are common 
to both PQ and unquenched correlators, the problem is reduced to counting contractions. 
The counting is the same for all planar quark-fiow diagrams with a given value of n, and 
thus eq. ()A1|) follows. 



APPENDIX B: CALCULATION USING ALTERNATE METHOD 

The aim of this appendix is to check the results, and to understand the limitations, of the 
primary method, by repeating the calculation using an alternate method, in which partial 
quenching is implemented using the replica trick. 

The starting point is eq. (jZj), which shows that one can obtain the average of the positive 
and negative spectral densities by determining the discontinuity in the twisted valence con- 
densate. To determine this condensate at LO using the replica trick one needs to minimize 
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the potential of a theory with Ns = 2 sea-quarks and 2Nv rephca quarks. Note that this 
theory is unquenched, and so has a well-defined and bounded Hamiltonian. The use of 2Nv 
replica quarks, rather than Ny, is for later convenience. Having calculated the conden- 
sate at the minimum, one then sends Ny 0, giving the valence condensate, and finally 
analytically continues to complex /i to determine the discontinuity across the imaginary axis. 

The LO chiral Lagrangian takes the same form as that in the PQ theory, eq. (pUj) . except 
that S G SU [Ns + 2Nv) and supertraces are replaced by traces. The potential is thus 



Tr + SxL'P) [Tr(S + St)]' [Tr(S-St)]' Tr (s^ + St 



Wq Wj Wg- 



P 4 16 16 16 

(Bl) 

The mass term takes the block form 

X±' = diag( 5 ^ i ^/"'''"S) • • • ,fn ± iJit^) . (B2) 

2x2 sea block A^^, 2x2 blocks 

The twisted mass appears only in the replica parts of the mass matrix since we are interested 
in the spectral density of Qm for a sea sector with an untwisted mass.^^ I recall that when 
using the rephca trick one takes the LECs {wq etc.) to be independent of Ny- 

I assume that the form of S which minimizes the potential follows that of the mass 
matrix, i.e. 

s = diag( ^ ,^^2^:^£2)- 

2x2 sea block AT^ 2x2 blocks 

This can be shown to lead to a local minimum under certain assumptions on the parameters, 
as I discuss at the end of the appendix. With the condensate given by ()B3|) the potential 
becomes 

^ ^„ = -m{l+Nvcosuj) - 'fiNvsinuj - Wg{1 + Nv cosoof - —{l + Ny cos2uj) (B4) 
/ 4 



constant — A''v 



rrf cos u; + u sin u; H cos 2u} 

4 



cos^ uj . (B5) 



Note that the w-j term does not contribute, and that the same combination rrf^ = m + 2wq 
appears as in the main text. As A''^ 0, the term quadratic in Ny can be dropped, and 



One can now see how the continuum resuh Q in the main text is obtained. Only the first term in V^°p 
is present, and x±' — diag(?7i5, my, . . . , fhv)- Minimizing the potential leads to the valence block of the 
condensate aligning in the direction of mv^ independent of the alignment in the sea-quark block. 
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the determination of the minimum is to be done using only the contribution hnear in Ny- 
This contribution is exactly the auxiliary potential encountered in the primary method, 
eq. (j56|) . Thus the required minimizations in the two methods are identical, and yield the 
same minimizing angle, ci;(/i). 

It remains to be shown that the results for the spectral density are also the same. In the 
replica method, the required condensate is, using eq. 



{qvlbT^qv) Lo = 2i5o/^ sin u . (B6) 

Since this is the condensate for one replica doublet (and not Ny doublets), it has a non- 
vanishing limit as Ny — > 0. Thus the final result from the replica trick is, using eq. (|7j), 

Bof Disc [smcu] = n [pq{\) + Pq(-A)] . (B7) 

li=iX 

Note that the spectral density pq is renormalized, since it is determined from a renormalized 
condensate. 

This result is equivalent to that from the primary method, eq. (|5n|) . To see this, one takes 
the derivative of the latter equation with respect to a, and uses the fact that the derivative 
commutes with taking the discontinuity: 

d 



-2mpR{a) = 4- { Disc [V^^^ipi'' = -z)] 



-Disc [^Vr (/i' = -^)] 



(B8) 



da 

Next, one notes that the spectral density of can be related to that of Qm, 

'l^/a pR{a) = pQ{^/a) + pQ{-^/a) , (B9) 



as well as the fact that 



= ^ = -2B,fsmu:{^^) . (BIO) 
c((/i^j dp, 



Substituting these results into eq. ()B8|) and performing some simple manipulations gives 
eq. ()B7|1 . The derivation also runs in the other direction: the result in the text can be 
obtained from ()B7|) as long as the discontinuity in sin a; is integrable, which is the case here. 

I now return to the question of whether the assumed form of the condensate, eq. ()B3|) . is 
correct. This brings up the general question of how to apply the replica trick in a context 
where the dependence on Ny is not perturbative. It seems possible, for example, that, 
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for a given choice of LECs wes, there could be a non-trivial Ny dependence of the phase 
structure, and thus a non-analytic dependence of the condensate on Ny- This issue deserves 
further study, which I have not attempted here. The only general comment I can make is 
that, if wq = Wj = 0, then, for wg > (which is the case in which the primary method fails), 
it is plausible that the condensate takes the form ()B3|) for all integer Ny and fi real. This 
is because the mass term wants to align S oc x+^j while the ws term is minimized when 
S = ±1. These competing "pulls" seem likely to lead to a condensate which lies on the 
shortest path between the two favored directions, which is eq. ()B3|) . Thus I do not expect 
that studying the phase structure will resolve the failures of the primary method in general. 

A more modest goal is to show that the form ()B3|) leads to a local minimum. It is 
straightforward to show that minimizing the restricted replica potential ()B5|) does lead to 
an extremum in all directions of possible fluctuations. Thus the issue is whether it is a 
minimum in all directions. Choosing a basis in which there is no mixing, I find, in the 
replica limit, that the pion mass-squareds are: 

Mls = m''' + ws, (Bll) 

Mhn = c (m^^ + ^J^) + ^8 {2c' - 1) , (B12) 

M^v,c = M^y^^ + wsil-c'), (B13) 

M'vs,n = M'yy^„ + wril-c)\ (B14) 

Mvs,c = ^ • (B15) 



Here c = cosu, the subscripts "c" and "n" stand for charged and neutral, respectively, and 
the other subscripts indicate the sea and valence content of the "pions". If any of these 
mass-squareds passes through zero and becomes negative, then the extremum is unstable 
to fluctuations in the corresponding directions. Note that the alternate method has more 
directions of fluctuation than appear in the primary method. Most notable are the "VS" 
directions — these introduce the LEG Wf which does not enter in the analysis of the main 
text. 

If the alternate method is to resolve the puzzle from the main text, it should do so for all 
choices of Wj, and so I begin with the simplest, Wi = 0. In this case, however, the results of 
the primary method are reproduced (as long as fh'^^ > 0, an exception discussed below). In 
particular, if wg > (and m'^^ > 0), all pion mass-squareds are positive for all real (positive 
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and negative) values of /i^, so that the choice ()B3|1 with u = u is a local minimum. If ws < 0, 
however, then Myy^ vanishes when /x^ = —a. This signals a non-analyticity in c, which is 
inherited by the valence condensate, and is identical to the condition for the gap discussed 
in the main text. Thus the puzzle is not resolved for these parameter choices. 

There are, however, two cases where the present analysis differs from that in the main text, 
both involving instabilities in which the mixed pions composed of valence and sea-quarks 
have vanishing masses. The most straightforward follows from eq. ()B14|) : as yU^ becomes 
more negative Myg ^ can pass through zero and become negative if is sufficiently negative. 
It turns out that should not be too negative, however, because then Myg^ also changes 
sign as /i^ becomes more positive. Nevertheless, the conclusion is that there is a range of 
negative which could lead to the expected analytic structure. 

The other possibility is more complicated. It occurs when rh^^ < 0, i.e when wq is 
sufficiently negative {2wq < —fh) that the mass in the effective potential is negative even 
though m itself is positive. Note that in this region one must have > \m'^^\, as shown 
in fig. ini As noted in the main text, if m''^ < 0, the minimum of auxiliary potential (or 
equivalently of V^'^'' for A'^v' 0) lies at c = — 1 if /i = 0. The mass formulae given above 
still hold, but now the contribution to Myg^ proportional to Wsc{l — c) is negative. This has 
the effect that Myg^ vanishes exactly at /i^ = 0, and becomes negative for /i^ < 0, for all 
allowed Ws- The vanishing at /x^ = can be understood from the viewpoint of the original 
replica potential, V^'^''. If c = —1, then the condensate takes the form 

(S) = diag(l, -!,...,-!). (B16) 

This breaks the SU{2 + 2Nv) vector symmetry of the theory that is exact when /i = 0. 
Thus, as in the Aoki-phase analysis, there are exact Goldstone bosons, here the VS pions. 
This symmetry breaking arises because, for wq negative, the wq term favors a condensate 
rotated away from ±1. 

This instability seems to imply a discontinuity in c, and thus in the valence condensate, 
beginning at — /i^ = 0. This would correspond to the gap vanishing. If correct, this could 
explain the failure of the primary method, which, unlike the present approach, relies on 
the gap being non-zero.^^ Note that one cannot reach the continuum limit in the region 

This illustrates a loophole in the claim made earlier in this Appendix that the two methods are equivalent. 
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of parameters where this possible resolution applies (m*^^ < 0, ws > \fh^^\), so there is no 
contradiction with the result that the gap is non-vanishing in the continuum. 

I have not pursued this analysis further, since the possible resolutions of the problem for 
Ws > are not general. In particular, they do not apply if m > > |2ty6|, which is the 
region through which one likely approaches the continuum limit if one is in the first order 
scenario. 
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